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Integral relaxation time of single-domain ferromagnetic particles 
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The integral relaxation time n n t of thermoactivating noninteracting single-domain ferromagnetic 
particles is calculated analytically in the geometry with a magnetic field H applied parallel to the 
easy axis. It is shown that the drastic deviation of r ; ~ t from the lowest eigenvalue of the Fokker- 
Planck equation Ai at low temperatures, starting from some critical value of H , is the consequence 
of the depletion of the upper potential well. In these conditions the integral relaxation time consists 
of two competing contributions corresponding to the overbarrier and intrawell relaxation processes. 
[S1063-651X(96)09409-3] 

PACS number(s): 05.40.+j 



I. INTRODUCTION 

At present, a single-domain ferromagnetic particle 
with uniaxial anisotropy attracts the attention of re- 
searchers, in particular as one of the models of informa- 
tion storage. The hysteretic rotation of the magnetiza- 
tion of such a particle over the potential barrier under the 
influence of an arbitrary directed magnetic field H was 
studied by Stoner and Wohlfart Jl]]. At nonzero tem- 
peratures the magnetization vector of the particle can 
surmount the barrier due to the thermal agitation, as 
argued by Neel this effect becomes especially pro- 
nounced for small particles having lower values of the 
potential barrier AU. Such a "superparamagnetic" be- 
havior was observed in many experiments on magnetic 
liquids, on polymers with magnetic inclusions, as well as 
on very thin magnetic layers forming "islands." 

An initial accurate calculation of the thermoactivation 
rate of a uniaxial ferromagnetic particle is due to Brown 
||, who derived the Fokker-Planck equation for an as- 
sembly of particles and solved it in the presence of a 
longitudinal magnetic field H = He Zl perturbatively in 
the low-barrier case, AU <C T, and with the use of the 
Kramers transition-state method [Q in the high-barrier 
limit T <C AU (the Boltzmann constant ks is set to 
unity). In both limiting cases considered by Brown the 
time dependence of the average magnetization (M z ) is a 
single exponential and the relaxation rate of ferromag- 
netic particles is given by the lowest eigenvalue Ai of 
the Sturm-Liouville equation associated with the Fokker- 
Planck equation. Subsequently, Ai was calculated nu- 
merically by Aharoni for arbitrary values of AU /T with- 
out a magnetic field j^j and with a longitudinal magnetic 
field |(| . The correction terms for the high-barrier result 
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for Ai were given by Brown Q. Later the analytical ex- 
pression for Ai in the high-barrier case was rederived in 
Ref. H with a more rigorous method. In Refs. §,|o|Jll) 
various approximate analytical formulas for Ai for the 
arbitrary AU /T were proposed. Recently the thermoac- 
tivation rate of single-domain magnetic particles, as de- 
scribed by Ai, was calculated numerically by Coffey et 
al. [[l2| for the arbitrarily directed magnetic field H, i.e., 
in the geometry considered by Stoner and Wohlfart (!]]. 

Apart from limiting cases, the Fokker-Planck equation 
for an assembly of single-domain ferromagnetic particles 
cannot be solved analytically. The magnetization relax- 
ation curve consists of an infinite number of exponentials 
and the overall deviation of the linear dynamic suscepti- 
bility from the Debye form can be as large as about 7% 
for isotropic particles in a static magnetic field, as shown 
in Ref. |13| . In this case it is convenient to introduce 
the so-called integral relaxation time T; nt , determined as 
the area under the relaxation curve after a sudden in- 
finitesimal change of the magnetic field. The quantity 
Ti n t depends on all eigenvalues A&, k = 1,2,..., and is 
therefore more informative than Ai; also it can be di- 
rectly measured. Moreover, it turned out that, unlike Ai, 
the integral relaxation time T; nt can be calculated analyt- 
ically for uniaxial particles in the longitudinal magnetic 
field for the arbitrary values of parameters |l3| and Tj"* 
recovers the analytical results of Brown for Ai in the 
asymptotic regions. 

The integral relaxation time was also the subject of 
a recent series of papers [PI JloTJlrl , where it was called 
the correlation time. In Ref. 14fl T; n t for uniaxial fer- 
romagnetic particles was calculated analytically with an 
alternative method for zero magnetic field, the resulting 
expression being, however, much more complicated than 
the original formula for T; n t of Ref. gjj). In Ref. 0] 
a numerical calculation of v m t in the case with nonzero 
longitudinal magnetic field was presented. In Ref. |l5| ] 
the congeneric model of rotating dipoles describing the 
dielectric relaxation was considered. The results of Ref. 
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|L4| show that in zero magnetic field tQ. is very close 
to Ai in the whole region of AU/T. On the contrary, 
numerical calculations of Ref. |16| reveal a striking be- 



havior r-~ t Ai for relatively small longitudinal fields 
in the region T <C AU. This region of parameters was 
not analyzed in Ref. |l3[ ], whereas in Ref. the effect 
was not physically interpreted. 

The aim of this paper is thus to consider in more de- 
tail the integral relaxation time of uniaxial ferromagnetic 
particles in the longitudinal magnetic field with the help 
of the method of Ref. |0|. As we shall see, the effect 
found in Ref. can be explained by the depletion of 
the upper biased potential well, which leads to the dom- 
inance of the fast relaxation inside the lower well in the 
integral relaxation time. 

The remainder of the paper is organized as follows. 
In Sec. [n] the derivation of the Fokker-Planck equation 
for an assembly of single-domain ferromagnetic particles 
from the stochastic Landau-Lifshitz equation is outlined. 
In Sec. Ill the known results for the thcrmoactivation 



rate of uniaxial ferromagnetic particles are briefly re- 
viewed. Then the integral relaxation time Tj n t is intro- 
duced and analyzed, and it is shown that the effect dis- 
covered in Ref. |L6| can be explained without an explicit 
calculation of r int . In Sec. [V the derivation of a gen- 
eral formula for r; nt of uniaxial particles in a longitudi- 
nal magnetic field is presented and its behavior is studied 
analytically and numerically for the whole region of pa- 
rameters. Some concluding remarks are given in Sec. 0. 



II. THE FOKKER-PLANCK EQUATION 

The magnetization of a single-domain ferromagnetic 
particle M can be considered not too close to the Curie 
point T c as a vector of fixed length |M| = M S (T), whose 
direction can fluctuate due to the thermal agitation. 
This fluctuative motion of M can be described semiphe- 
nomenologically with the help of the stochastic equation 



M = 7 [Mx (H off + C)] -R(M), 
where 7 is the gyromagnetic ratio, 



Hcff -~aM' 



w 



H M - km: 



(2.1) 



(2.2) 



are the effective field and the energy density, H is the 
external magnetic field, and K is the uniaxial anisotropy 
constant. The energy of a particle is given by 



U = VW, 



(2.3) 



where V is the particle volume. The correlators of dif- 
ferent components of the white-noise field £(i) can be 
conveniently written as 



(Ci(*)0(*')> 
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(2.4) 



The relaxation term R in (2.1) describes, like C, the in- 
fluence of the heat bath on the particle and, as we shall 
see immediately, it has the Landau-Lifshitz form |l7j 



R = 7 A[M x [M x H cff ] 



(2.5) 



The Fokker-Planck equation corresponding to ( |2.1| ) is 
formulated for the distribution function /(N, t) = (i5(N— 
M(i))) on the sphere |N| = M s , where the average is 
taken over the r eali zations of £. Differentiating / over t 
with the use of (2^) and calculating the right-hand side 
analogously to the derivations given, e.g., in Refs. Jl8| , p^| , 
one comes to the Fokker-Planck equation 
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(2.6) 



One can easily see that the equilibrium distribution func- 
tion 



/o(N) oc exp[-W(N)/T], 



(2.7) 



is the solution of fl2.6|) if and only if R has the double- 
vector product form ( |2.5| ) , which reflects the way of how 
magnetization is coupled to the heat-bath fluctuations in 
fl2.ip . If, e.g., the correlators of C components i n (|2.4 |) 
are anisotropic, the expression ( |2.5| ) also changes |13||. 

Brown used in his derivation of the Fokk er-P lanck 
equation || the stochastic equation of motion ( |2.l| ) with 
the Gilbert relaxation term R = 777 [M x M] pffl. Re- 
defining 7 =>• 7g in Brown's equation, one can trans- 
form the latter to the form (2.1) and (2.5) with 7 = 



7g/(1 + V ig^s) an d ^ = VJg, where the Langevin 



field C enters also the expression for R (2.5) as being 
added to H c fj. This means using a stochastic model 
somewhat different from the one described above. Both 
models coinside, however, in the actual small-damping 
case rf)M s <C 1. 

The equation of motion for the magnetization (M) = 
J <i 3 A^N/(N, t) o f an assembly of particles can be easily 
derived from (|2.6|) and has the form 



> 



7 ([M x H off ]) 
-A([M x [M x H cff ]]) - A N (M) 



(2.8) 



[cf. (2.1)], where the characteristic diffusional relaxation 
rate Ajv is given by 



A 



•N 



2-yXT/V. 



(2.9) 



One can see that even in the case without anisotropy, 
where H e g = H, this equation is not closed since it is con- 
nected to the second-order corr elati on functions (MiMj) 
in the Landau-Lifshitz term of (|2.8| ). Therefore, the res- 
onance a nd r elaxational behavior of the Fokker-Planck 
equation ( |2.6| ) is in general not described by Lorentz and 
Debye curves, and the deviations from the latter can be 
about 7% Q. Neglecting these features, one can ob- 
tain the best isolated equation of motion for the magne- 
tization of an assembly of particles in the isotropic case 
K = 0, choosing the distribution function in the form 



2 



/(N,i) ocexp[yA(i)N/T] [cf. (|2.7|)1, where the temporal where h is given by (^3|). The factor (1 + h) before 



Such a 
contains 



evolution of the vector A(£) is governed by (2 
generalized Landau-Lifshitz-Bloch equation 
both transverse and longitudinal relaxation terms. In 
the high-temperature limit (in the isotropic case K = 
this requires T ^> VHM S ) Eq. ( |2.8| ) becomes closed and 
takes on the form of the Bloch equation with the relax- 
ation rate Ajv- 



III. THE LOWEST EIGENVALUE Ai AND THE 
INTEGRAL RELAXATION TIME n NT 

To parametrize effects of thermal agitation on ferro- 
magnetic particles, it is convenient to introduce the di- 
mensionless energy u = Ji/T = VW/T, which in the case 
with a longitudinal magnetic field has the form 



x = cosd = M z /M s 



with 



. VHM S VKMg 

£ = — ^ — i a — — — • 

The top of the barrier corresponds to 
x - — - x - — - ■ h . h — ~~ 



(3.1) 



(3.2) 



(3.3) 
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The barrier height Am = u(x m ) — u{— 1) is given by 

Au = a-Z + £ 2 /(4a) = a(l-h) 2 . (3.4) 

In the case a,£ ~ 1 a general solution of the Fokker- 
Planck equation ( |2.6| ) cannot be found analytically and 
the relaxation of any initial state is described by a sum 
of exponentials of the type A exp(— Aji), where A, are 
the eigenvalues of the Sturm-Liouville equation associ- 
ated with the Fokker-Planck equation (for the longitu- 
dinal relaxation all A^ are real). In the low-barrier case 
a,£ <C 1, the problem can be solved perturbatively ||] 
and the longitudinal relaxation is governed with a good 
accuracy by the single exponential corresponding to the 
lowest eigenvalue Ai which is given by 0] 



Ai = A 



2 

1 — -a 
5 
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(3.5) 



with A^v determined by ( |2. 

In the high-barrier case q>1, the relaxation is dom- 
inated again by Ai describing now the slow overbarricr 
thermoactivation, whereas all other eigenvalues Aj cor- 
respond to the fast intrawell processes with small am- 
plitudes. Brown's result for the high-barrier case, which 
was derived with the help of the transition-state method 
of Kramers [Q , can be written in the form || 

Ax ~ ANir-^a^il - h 2 ) (3.6) 
X {(1 + h) exp[-a(l + h) 2 ] + (1 - h) exp[-a(l - h) 2 ] }, 



the exponenti al fu nction in (3.6) is irrelevant since the 
first term of (3.6) is only essential for £< 1, which for 
a ^> 1 implies h <C 1. We will, however, keep this factor 
here and in analogous expressions below for the sake of 
symmetry. 

In the intermediate region a, £ ~ 1, it is convenient 
to introduce the integral relaxation time determined as 
the area under the magnetization relaxation curve after 
a sudden infinitesimal change of the applied field H by 
AH at t = 0: 



Tint 



(It 



(M z (oo)) - {M z (t)) 
(M Z (^))-(M Z (Q))- 



(3.7) 



Unlike Ai , the integral relaxation time Tint ca n be found 
analytically from the Fokker-Planck equation (|2.6|) in the 
whole range of parameters in the geometry with a longi- 
tudinal magnetic field Jl3| , as will be described in detail 
in Sec. IV. Here we discuss the results of recent calcu- 
lations of Tint by Coffey et al. jl4],[l6]]. At first note that 
the relaxation curve can be represented in the form 



(M,(oo)) - (M z {t)) = AHxz^Ai 



-Ait 



(3.8) 



where \z = d(M z ) / dH is the static longitudinal suscep- 
tibility. This form of writing the response function is 
more convenient than that of Refs. |l4|, 16| since here the 
amplitudes Aj obey the sum rule Ai = 1. Now t 1vA of 
J\j can be rewritten as (cf. Ill]) 



Tint 



(3.9) 



In Refs. |14|,[ll| the integral relaxation time Tj nt is called 
the correlation time since according to the fluctuation- 
dissipation theorem r; n t can be also considered as the 
area under the autocorrelation function. The term "cor- 
relation time," however, seems to be rather artificial be- 
cause the autocorrelation function does not appear in the 
actual calculation of T; nt with the help of (3.7) or ([3.9]), 



as well as in Sec. IV below, and really considering auto- 
correlations would imply going unnecessarily beyond the 
Fokker-Planck equation. 

According to the numerical results of Coffey et al. 
|l4j in a zero magnetic field the amplitudes A; satisfy 
Ai -C A\, i = 2, 3, . . ., for all values of a and the differ- 
ence between Ai and r^ 1 is small everywhere reaching 
only 1.2% at a = 5. On the contrary, the subsequent 
calculations for H =/= revealed a striking behavior 
T int A-i at sufficiently low temperatures. The formal 
reason for this is that A\ becomes small in this region 
and the terms with k — 4, 5 dominate in ( |3.9| ), as shown 
in Ref. Jjl|. But the effect can also be interpreted on a 
physical level as the consequence of the depletion of the 
upper potential well and quantitatively described with- 
out a general calculation of T; nt , as will be demonstrated 
below. 

The reduced equilibrium magnetization of an ensemble 
of noninteracting ferromagnetic particles m z = (M z )/M 3 
is given by the generalized Langevin function -B(£, a): 
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f 1 d 
J xf (x)dx = —\nZ = B(£_,a), 



(3.10) 



where, according to (2.7) and (3.1) 



/o = ^exp(-W/T) = ^ 



Z= I e- u dx. (3.11) 



In the high-barrier case a>l, the partition function Z 
is a sum of two contributions corresponding to the two 
potential wells Z = Z + + Z_, 



Z ± ^ 



2a ±£ 



2a 



(2a ±0' 



(3.12) 



where the correction terms account for the curvature of 



the potential-energy function u of (3.1). Neglect ing th ese 
small terms, one can represent B(£,a) of Eq. (3.10) by 
two mutually complementing expressions 



B(£,a) =tanh£- — ' 



2a Vcosh 2 ^ 



tanh£ 



for 1 ~ £ <C a and 



B{£,a) = 1 - 2e 



^^±z ]_ 



2a - £ 2a + £ 



(2a)* 
(3.13) 



(3.14) 



for 1 <C £ ~ a. Here, in the first limiting expression 
the second term is small and irrelevant; the third term 
is kept since it yields a contribution to the derivative 
B' = dB/dt; that is not exponentially small for £ 3> 1. In 
the second limiting expression (the strong-bias case) the 
deviation of B from unity separates into two parts: The 
second exponentially small term is due to the population 
of the upper well (x ~ — 1), whereas the third one ac- 
counts for the thermal agitation in the lower well (% ~ 1). 
The response of magnetization to an infinitesimal change 
of the magnetic field AH is related to the derivative of 
the generalized Langevin function B': A m = B'AH /T. 



The latter can be determined from (3.13) and (3.14) and 
put in the whole region 2a — £ ^> 1 into the unique ex- 
pression 



B' = B'p, + B' w = t\ — t 



1 



cosh 2 £2a-£ ' (2a + 2 
1-h 2 1 



ci&h) (2a + 2 



where 



^h) = -[(l + h)e-' + (l-h)e^ 



(3.15) 



(3.16) 



and B' B accounts for the redistribution of particles be- 
tween the two wells across the potential barrier and B' w 
that inside the lower well. Henceforth we will use the sec- 



ond of the equivalent forms of B' in (3.15) for the sake 
of symmetry [cf. the comments after Eq. (3.6)]. 



mJ°°)-m z (f), arb. units 



The magnetization relaxation curve 
in the high-barrier strong-biase case: 
the intrawell (1) and overbarrier (2) 
contributions, 
Am B «Am w , A l «A w 




1 - Am w exp(-?A w ) 

2 - Am B exp(-;Aj) 

3 - the total relaxation curve 



6 8 10 

t, arb. units 

FIG. 1. A schematic look of the two-exponential relax- 
ation curve of single-domain ferromagnetic particles in the 
strong-biase high-barrier case, «, £ S> 1. 



Now, in the low-te mper ature strong-bias case a, £ 3> 1, 
the relaxation curve (3.S) consists of only two exponen- 
tials (see Fig. |l|) 



m 2 (oo) - m z (t) 

= AniB exp(— tAi) + Amw exp(— tAw), (3-17) 



(3.18) 



where Ai is given by ( |3.6| ) and 

A w = 2 7 Aff c ff = A N (2a + f) 



is the temperature-independent relaxation rate in the 
lower well, which can be obtained from the determin- 



istic Landau-Lifshitz equation (2.1) and (2.5) without £. 
The integral r elax ation time Tint calculated according to 
the definition (3.7) can be written as 



Tint = Tint R + U n t W, 



where 



Ti n t,B 



B 



B a -1 



Tint,W = 



B' 



W A-l 



(3.19) 



(3.20) 



One can see that in the low-temperature strong-bias case 
a,£ ^> 1, the barrier contribution Ti n t,s into the integral 
relaxation time Tint can be substantially reduced due to 
the depletion of the upper potential well manifesting itself 
in the exponential smallness of the magnetization change 
due to overbarrier transitions Arris oc B' B [see (3.15) 



In this case the overbarrier and intrawell terms in (3.1£) 



can compete with each other since Ai is exponentially 
small and B' w /B' = 1. On the contrary, for small or 
zero bias one has B' B /B' = 1 and B' w /B' <C 1, so that 
the intrawell process can be completely ig nored . 

The expressions for T; nti s and T m t,w in (|3.20| ) are valid 
in the whole high-barrier region 2a — £ 3> 1 and will 
be obtained independently in the framework of a general 
method in Sec. IV. In the strong-bias case £ 3> 1, the 
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barrier contribution rj n t,B in ( |3.2C|) c an be represented 
with the use of (pQj), |tl~5| ), and (|3.3[) as 



IGA^^q) 1 / 2 ^^ exp[a(l -6h + h 2 )). 



(3.21) 



It changes its behavior as a function of a at the critical 
value of the applied field 



h = h c = 3- 2\/2 w 0.17, 



(3.22) 



which is substantially smaller t han the fiel d o f the bar- 
rier disappearance h — 1 [see (3.3) and (3.4 )1. F or h 
in the vicinity of h c the exponential factor in ( 3.21 ) can 
be written as exp[— A^/2a(h — h c )\. It can be seen that 
for h < h c the quantity Ti n t,B exponentially increases 
with lowering temperature (i.e., with increasing of a) an d 
brings the dominant contribution into T m t of Eq. (3.20). 
On the contrary, for h > h c the quantity Ti n t,B expo- 
nentially decreases at large a, so that Tint tends to the 
temperature-independent value T m t,w of (3.18). One can 



also see that for h only slightly higher than h c the quan- 
tity T- lnt ,B increases as a 1 / 2 at smaller a; then the de- 
creasing exponential becomes to dominate. Thus, in this 
case Tint^s, and hence Tint of (3.20), has a maximum at 



some a»l and r ; ~ t has the corresponding minimum, as 
was obtained numerically in Ref . jl6| . It should be noted, 
however, that the actual position of this min imum can be 
described only takin g into account in ( 3.20| ) the general 
form of B' given by ( [3.15 ) . 

The results above completely describe the observations 
made in Ref. |l6[ ] in the low-temperature strong-bias re- 
gion. In the next section we present the analytical calcu- 
lation of ri n t in the whole range of parameters. 



IV. CALCULATION OF THE INTEGRAL 
RELAXATION TIME tint 

All the information about the relaxation curve ( |3.8| ) 
is contained in the longitudinal linear dynamic suscepti- 
bility x(w). In the presence of a small alternating field 
AH z (t) = AH z q exp(— iujt) the deviation of the distribu- 
tion function / from the equilibrium function (3.11) can 
be represented as 



Sf = f (x)q(x)VM s AH z (t)/T, 



(4.1) 



where the function q(x) satisfies an equation following 
from (|27 



— + 2ax 

ax 



2 dq 



2iwA 



n 1 



= (l-x 2 )(2ax + £)-2x. 



(4.2) 



The dynamic susceptibility of the particle's assembly is 
then determined by 



VM* 
T 



dxxf {x)q(x). 



(4.3) 



Using the linear-response theory one can easily show 
that Xz( w ) can be represented in the form 



A, 



(4.4) 



where Xz-, A i7 and A 4 are t he parameters of the mag- 
netization relaxation curve (3J3). Calculating this sum 
of Debye terms requires knowing all eigenvalues Aj and 
amplitudes Ai associated with the Fokker-Planck equa- 
tion and cannot b e do ne analytically in the general case. 
Accordingly, Eq. ( |4.2| ) has no general analytical solution 
and its behavior is to be studied analytically in the lim- 
iting cases of high and low temperatures and high and 
low frequencies, as was done in Ref. p3| . In particular, 
generating high-frequency expansions of x( w ) does not 
require solving differrential equations and can be carried 
out up to high orders. The corresponding results, how- 
ever, are not very interesting here since they describe 
only fast intrawell processes. The information about the 
slow process of thermoactivation is contained in the low- 
frequency expansion of x(w), which can be written in the 
form 



(4.5) 



with ( |4.4| ) , one can show that the quan- 
tity ri n t in (4.5) is exactly the integral relaxation time 



Comparing (4.5 



given by the formula ( |3.9| ) 

The perturbative solution of (4.2) for small u> can be 
done analytically since for to = there are only terms 
of the type q'(x) and q"(x) in the equation. Hence one 
can introduce a new variable g{x) = q'(x) and solve suc- 
cessively the first-order differential equation s for g(x) and 
q(x). After calculation of the susceptibility (4.3) one gets 
the analytic expression the for integral relaxation time 

Tint 0: 



Tint 



dx 



AjvB' J_ x 1 - x- 



where f is given by (|3.1l| ), B' — dB/d^, and 



(B - x')f (x')dx' 



(4.6) 



(4.7) 



Recalling the general formula for B(£, a) Eq. ( p.lOD , one 
can conclude that $(±1) = 0, i.e., the integrand of (4.6) 
goes to zero at x = ±1. The function <i>(x) can be easily 
calculated analytically in two particular cases. In the 
unbiased case £ = one gets 

= 7?-[/o(l) - fo(x)], fo(x) = e ^p. , (4.8) 
2a Z(a) 



whereas in the isotropic case a = 



$(x) 



fo(x) 



coth £ — x 



exp(-£z) 
sinh £ 



(4.9) 
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with fo(x) — £exp(£x)/(2sinh£). An alternative analyt- 
ical expression for T m t in the particular case £ = was 
derived recently by Coffey et al. jl4j with the help of the 
development of the solution of the Fokker-Planck equa- 
tion in Legendre polynomials. Their expression contains 
Kum mer functio ns and is essentially more complicated 
than (fD| with (p~g|). 

Now w e pr oceed with the analysis of the general ex- 
pression (4.6) in different limiting cases. In the high- 
temperature limit a, £ <C 1, the calculation can be done 
perturbatively with respect to a and £. In particular, 
at very large temperatures one has fo(x) = 1/2, B = 0, 
B' = 1/3, and $ = (1 - x 2 )/4. Thus the whole phase 
space of the fer romagnetic particle, — 1 < x < 1, con- 
tributes to (4.6) and one gets T;",. 1 = Ajy. A more accu- 
rate calculation yields 



15 



10 



Zd> 2 (x)/[f Q (x)(l-x 2 )] 



20 h = 0.35 

a = 2a+t, = 


0, 5, 10, 15, 20 
20 


/ 15 \ 




/ / 10 \\ 




//^--~~-~~~~~~~~ o 


/ 15 



-0.5 



0.0 



0.5 



A 



N 



1 



35 



(4.10) 



whic h is very close to Brown's expression for Ai given by 

In the unbiased low-te mper ature case £ = 0, a 1, 
the function $ given by (4.S) is constant in the main 
part of the x interval, except for near t he b orders. Thus 
the main contribution to the integral (4.6) comes from 
the barrier region x ~ a" 1 / 2 -C 1 cut by the function 
Jq 1 (x exp(— ax 2 ). With the use of (3.12) with £ = one 
gets, in the leading order, 



- 2A A ^- 1 /2 a 3/2 £ 



which coincides with the expression for Ai in (|3.6| ) in 
the unbiased case h = 0. It is not difficult to calcu- 
late also the correction terms for the formula ( 4.11 ), 
which coincide with those given by Brown |Q. In the 
isotropic strong-field limit a = 0, £ ^> 1, the function $ 



given by (IE), as well as the whole integrand of (4.6), is 
peaked in the vicinity of the potential minimum x = 1 , 
where the exponentially small term with exp(— £x) can 
be neglected. In the leading order o ne ge ts for tQ the 
temperature- independent expression ( 3.18|) with a = 0. 

Now we consider, as a corollary, the biased low- 
temperature case $ / 0, a > 1. As argued in Sec. ID, 
in this case Tj nt can be the sum of intraw ell an d overbar- 
rier contributions Ti n t t w and [see ( 3.1E)| )1. Accord- 

ingly, the integrand of fl4.6| ) can consist, for a, £ 3> 1, 
of two peaks (see Fig. 2|) corresponding to the barrier 
top x ~ x m — —h in ( 3.3| ) and the lower well x *~ 1. 
The function $(x) of ( [4.7D is determined by two well- 
separated potential wells and is therefore practically in- 
dependent of x for x not too close to the boundaries. In 
the the lower- well region 1 — x <C 1, the calculation yields 
$(x) as a sum of two contributions <& = $b + $vy , where 
to the leading order 



$ B (x) 



1-/i 2 

2c 2 (e,M 



{l-exp[-(2a + £)(l-x)]}, (4.12) 



c(£, /i) is given by (3.16) and 



FIG. 2. The integrand of the formula ( |4.6|) for the integral 
relaxation time Tint for different values of parameters. 



$ w (x) =(l-x) exp[-(2a + £)(1 - x)]. 



(4.13) 



The term $b of ( 4.12| ) goes over to the constant men- 
tioned above in the region not too close to the border 
[1-x > l/(2a+£) -C 1]. For £ > 1 it acquires the small 
factor e -2 ^ acco untin g for t he de pletion of the upper po- 
tential well [cf. ( 3.15 ) and ( 3.20 )]. On th e contrary, such 
a factor is not present in $vy m (4.13), but the corre- 
sponding contribution into T; nt given by (4.6) is reduced 



(4.11) due to / c 



gets Eq. ( 3. IE ) , where 



Now calculating the integral (4.6) one 



B' 



rV2 



Tint.B 



exp[a + £ 2 /(4a)] 



B> 2A N a 3 / 2 (l-h 2 )c{£,h) 



(4.14) 



coincides with the expression given by ( |3.20 ) and in the 
strong-bias case 



T mt,W 



A 



w 



An 2a + £ - 



10a + £ 
2a + £ 



(4.15) 



In ( 4.15| ) [cf. ( 3.18 )1 the c orrection terms, in particular of 
the type present in ( |3.12] ), have been taken into account. 

Numerical c alcu lation of the integral relaxation time 
Tint given by (4.6) in the whole range of parameters a 
and £ poses no difficulties. For the representation of the 
results for the arbitrary relation between a and £, in- 
cluding the case a = 0, it is more convenient to use the 
variable 



a = 2a + £ = — —wr, 



7 T 



(4.16) 



where ujr — r )(2KM s +H) is the ferromagnetic resonance 
frequency in the lower potential well. One can see that 
2a = a/(l + h) and £ = ah/ (I + h). In term s of the 
variables a and h the asymptotic formula ( 4.15] ) can be 
rewritten as t^ w = Ajv [a — (5 + h)/ (1 + h)] , which 
shows that all curves t^ w (cl) for different h are parallel 
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FIG. 3. Inverse of the integ ral r elaxation time Tint as a 
function of the parameter a of (4.16). Circles: data of Coffey 
et al. jl(| taken from their Tab. II. 



to each othe r. T he results of the numerical integration in 
the formula ( |4.6| ) are represented in Fig. |[ These results 
confirm all the considerations made above, as well as the 
numerical findings of Coffey et al. |Iq ]. 



respond to the intrawell and overbarrier processes. In 
this case it would be naturally an oversimplification to 
describe the relaxation with a common integral relax- 
ation time, but it is possible to separa te a nalyti cally both 
contr ibutions in the general formula (4.6) [see ( 4.14 ) and 
(4.15)] and thus to give a complete description of the 



relaxation process. 

Going beyond the integral relaxation time, one can 
conceive the calculation of the wh ole dynamic suscep- 
tibility Xz(w) given by the formula (4.4). In this case the 
numerical approach of Ref . Jl6| , consisting of the calcu- 
lation of all Ai and Aj, is really useful. However, for the 
model with a high potential barrier such a calculation is 
probably not very interesting for the reasons mentioned 
above: there are, to a good accuracy, only one (in the 
unbiased case) or two (in the strong-bias case) terms in 
(4.4). More appealing would be to produce calculations 
for an isotropic model in a field (a = 0, £ ^ 0) where the 
eigenvalues Aj are not so well separated. In this case de- 
viations from the Debye form of the longitudinal dynamic 
susceptibility were studied analytically in Ref. jl3) , where 
it was shown that these deviations can be about 7% at 
(«3. In Ref. jl3| an effective two-relaxator formula for 
Xz(^) was proposed, which was argued to describe the 
main part of deviations from the simple Debye form. It 
would be interesting to check and improve these results 
by a direct numerical calculation. 

One more unsolved problem is the calculation of the 
transverse integral relaxation time of superparamagnetic 
particles. Up to now it seems to be considered only for 
the model of rotating dipoles in Ref. ]l5| . 

Although many experimental investigations are cur- 
rently done on systems showing superparamagnetism, 
these investigations are practically confined to the cer- 
tification of a superparamagnetic behavior and to rough 
estimation of relaxation times. It would be worth making 
more purposeful measurements aimed at a comparison 
with existing theories. For this purpose it would be im- 
portant to eliminate the distribution of particle volumes 
and orientations of the anisotropy axes. 
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tential well the relaxation is described by two exponen- 
tials with uncomparable characteristic times, which cor- 
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